Spin filtering by ferromagnetic nanowires 
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We show that electrical current flowing through nanowires made of ferromagnetic disordered alloys 
can become highly spin polarized. 
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The increasing ability to control spin-dependent phe- 
nomena in condensed matter physics is opening exciting 
possibilities for the electronic industry. The search for 
efficient spin filtering systems in particular is very im- 
portant for spin-based electronics and quantum compu- 
tation. Spin polarized electrical current for example may 
be required to operate qubits in some proposed device 
schemes for quantum computations jl], g, [| . 

Recently, a clever and promising technique has been 
devised for fabricating very thin metallic nanowires, with 
diameters smaller than lOnm H. It uses suspended car- 
bon nanotubes as substrates for deposition by electron- 
beam evaporation of several metals, including Fe and Ni. 
By initially coating the carbon nanotube with a small 
amount of Ti (with nominal thicknesses ranging from 1- 
2nm) other metals can stick more easily to the substrate, 
leading to the formation of nice continuous nanowires up 
to tens of microns long. 

The purpose of this letter is to show that electrical 
current flowing through nanowires of ferromagnetic dis- 
ordered alloys can become highly spin polarized. There- 
fore, wires with such characteristics may be useful as spin 
injectors in electronic devices. We begin by presenting 
and discussing the main physical mechanisms that are re- 
sponsible for the spin filtering behavior of these systems. 
Then we perform model calculations to illustrate how ef- 
fective those mechanisms can be, and finally discuss our 
results and the guidelines they provide for choosing the 
nanowire composition that may maximize the spin filter- 
ing effect. 

It is well known that disorder in metallic systems may 
lead to localization. Depending on the nature and de- 
gree of disorder, localized states may appear near the 
top and bottom of the conduction bands of three dimen- 
sional metallic systems, with a mobility edge separating 
localized from extended states. Non-interacting disor- 
dered electronic systems with lower dimensionality in the 
thermodynamic limit have localized states only. Strictly 
speaking they are insulators at zero temperature [j5], [|. 
In fact, localization effects can be very effective in reduc- 
ing the conductance of nanowires ffij. They can make 
the average conductance g fall off exponentially with the 
wire length £ for sufficiently long wires. The localization 
length A is determined by the asymptotic decaying rate 



of g, being shorter the faster g decreases. Nikolic and 
MacKinnon have made a detailed study of the electrical 
conductance in non-magnetic disordered nanowires, em- 
ploying a single-band tight-binding model to describe the 
electronic states |tJ . The problem involves basically four 
characteristic lengths: £, A, the wire width w, and the 
electronic mean free path A. In their work, they illus- 
trated the occurrence of different transport regimes ac- 
cording to existing relations between £, A, and A. Trans- 
port is quasi ballistic, when I is comparable with A, meso- 
scopic when A < £ < A, and strongly localized when 
A < £. The quantities A, A, and g are all functions of 
energy, and the conductance may change regime as the 
energy varies. Actually, in the presence of bulk disorder 
they found that g decays faster as a function of £ for en- 
ergies close to the band edges in the strong localization 
regime. 

It is noteworthy that the scattering of carriers in metal- 
lic ferromagnets is generally spin dependent, even when 
the scattering potentials do not depend upon spin. This 
is mainly due to the densities of states around the Fermi 
energy being different for majority and minority spin car- 
riers. In ferromagnetic transition metals, for example, 
the spin polarization of the sp-electrons is relatively small 
compared with that of the d-electrons. Both sp and d- 
electrons participate in electrical conductance, but the 
d-electrons are less mobile because they have a larger ef- 
fective mass. Nevertheless, even assuming the current is 
predominantly carried by sp-electrons in such systems, 
sp-d hybridization may lead to distinct conductances for 
up and down spins. The reason, as rightfully argued by 
Mott, is that electrical resistance is proportional not only 
to the density of scattering centers, but also to the num- 
ber of available states where electrons can scatter into. 
Thus, the existence of unoccupied d states at the Fermi 
energy [Ep) in transition metals acts as a trap for the 
sp electrons, because sp-d hybridization allows them to 
be scattered into the available d-states. Since the den- 
sities of available d-states at Ep differ for up and down 
spins in ferromagnetic transition metals, it follows that 
the electronic mean free paths of majority and minority 
spin carriers are usually not equal in such systems. Fur- 
thermore, the atomic potential fluctuations experienced 
by d electrons in transition metal alloys are often much 
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larger than those felt by s electrons. As a consequence, s 
and d electrons are differently affected by disorder, and 
we will show that the localization lengths for majority 
and minority spin electrons may be also rather different 
in ferromagnetic nanowires made of transition metal al- 
loys. In such systems A, and A, are both energy and spin 
dependent quantities. 

An additional relevant length scale for discussing spin 
dependent transport in ferromagnets is the mean free 
path associated with spin flip scattering (Xj). In metallic 
systems A/ is usually much longer than A. For nanowires 
with £ < Xf , the conductances for up and down spins are 
independent, hence g = g~] + g~i, where g a is the aver- 
age conductance for electrons with spin a. The up- and 
down-spin channels in this case behave as resistors in par- 
allel. Since both g^ and depend upon energy, they may 
decay as functions of l with different rates, particularly 
in strong ferromagnets where the majority d-bands are 
completely filled. With dissimilar localization lengths for 
up and down spin carriers, the polarization of the elec- 
trical current may thus increase very rapidly with t. 

To illustrate how effective these mechanisms can be, we 
perform model calculations for ferromagnetic nanowires 
made of metallic disordered alloys. For numerical 
simplicity we consider two-dimensional wires of finite 
lengths, sandwiched by perfect leads. Experimentally the 
wire thickness is not constant, and this is simulated by 
assuming it fluctuates randomly by ±5w around an aver- 
age width w. In order to allow spin dependent Mott scat- 
tering mechanism to take place, the nanowire electronic 
structure is described by a simple s-d model Hamiltonian 

where aj is the creation operator for an electron with 
spin a in orbital n on site i, and is the corresponding 
occupation number. The sp- and d-orbital sets in tran- 
sition metals are represented here by just two s orbitals 
labeled by ji = 1, 2, respectively. The on-site matrix el- 
ements h?" ~ tiSftu + (1 — &nv)li, where are atomic 
energies, and ji characterizes the hybridization between 
the s and d bands. Our system consists of a ferromag- 
netic disordered alloy, where sites are occupied with a 
certain probability either by magnetic or non-magnetic 
metal atoms. We neglect the effect of disorder in the 
hopping integrals, and consider them as being non zero 
between nearest neighbor sites only. Thus, for nearest 
neighbor sites i ^ j, hy" — —t^S^, where t ss and t dd 
symbolize the transfer integrals for s and d electrons, re- 
spectively. U- 1 represents the Coulomb interaction be- 
tween electrons located on the same site and orbital, and 
we further assume it takes place only when they occupy 
the d-orbitals of the magnetic metal atoms. 

In principle all the parameters involved here can be es- 
timated from the band structures or renormalized atomic 



potentials of the constituent metals. However, given the 
model nature of our calculation we have selected repre- 
sentative values for the parameters, rather than attempt- 
ing to adjust them to fit a specific system. For example, 
we choose our energy unit such that t ss — 1, and take 
t dd = 0.2, and ji = 0.9 independent of the lattice sites. 
This approximately reproduces the typical bandwidths 
of transition metals. By placing the band structures of 
the constituents on the same absolute scale, and aligning 
the corresponding Fermi energies, one may also estimate 
the energy levels characterizing the host (H) and 

impurity (I) metals. We shall discuss alloys based on 
ferromagnetic transition metals such as Fe, Co or Ni, 
with non-magnetic transition metal impurities of the left 
side of those elements in the periodic table. We thus 
set e s H = as our energy origin, and choose e d H = 0, 
ej = 0.03, and e d — 0.6 or 0.9, as a reasonable set of 
values for describing representative systems in the scope 
of such simple a model. 

We treat the electron interaction within the Hartree- 
Fock approximation, thus reducing the on-site interaction 
term to — hA^d^a, where represent the exchange 
splittings for the s and d orbitals, and a = ±1 for f and 
J, spins, respectively. With the simple form of interaction 
assumed, only A d H ^ 0, corresponding to a molecular 
field acting solely on the d orbitals of the magnetic-metal 
host atoms. 

Since our wire may be viewed as a sequence of atomic 
chains of finite sizes w ± 5w, it is convenient to label the 
atomic site positions by a pair of indices (Z,n) represent- 
ing the line / the atom belongs to, and its position r\ 
along that line. The conductance in the spin channel a 
is calculated by the Kubo formula || 

4e 2 

9a{E) = —— ReTrfGgo^oiG^tio — toiG1 Q toiG1 Q ] , (2) 

where and 1 symbolize any two adjacent line indices 
(the choice is arbitrary due to current conservation), and 
tax is the tight-binding hopping matrix between such 
lines. Gf m = &(Gf- - Gf+), where G\% are matri- 
ces representing the advanced and retarded one-electron 
propagators for particles with spin tr, evaluated at en- 
ergy E, connecting lines I and m. ReTr stands for the 
real part of the trace over all orbitals and line sites, and 
e /h is the quantum conductance. 

The required one-electron Green functions were calcu- 
lated by the method described in reference ||. Firstly, 
the surface Green functions of the semi-infinite perfect 
leads are generated by a well established technique. j|] 
Then, the Dyson equation is employed recursively to built 
the corrugated disordered wire parts atop the leads to the 
left and right of lines and 1, respectively. Finally, those 
wire parts are reconnected by turning on the electron 
hopping toi between them. In our calculations the impu- 
rity sites are randomly chosen for a given alloy concen- 
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tration, and configurational averages of g a are performed 
over 4,000 different samples. 

Before presenting our results for ferromagnetic sys- 
tems, it is instructive to study the conductance of non- 
magnetic disordered wires with this simple s-d model, 
particularly the behavior of g as a function of I for differ- 
ent values of energy. Thus, we start by setting A d H = 0, 
and consider wires with w — 10 atoms, Sw = ±1 atom, 
sandwiched by leads of width wl — w (leads of larger 
widths play no significant role on the main points we 
wish to address). Figure 1 shows localization lengths cal- 
culated as functions of energy for different impurity po- 
tentials and concentrations. Here, it is statistically more 
appropriate to define g = exp[(lng}], where (/) denotes 
the configurational average of / J/J. A is determined by 
fitting the asymptotic exponential decaying rate of g with 
i; more precisely, A -1 = —dg/dl for I ^> 1. The localiza- 
tion length clearly depends upon the impurity potential 
strength relative to the host, i.e., on Ae d = ef — t d H . It 
is also evident that A -1 increases as energy approaches 
the d-bands region (delimited here by « e d ±0.8 for non- 
hybridized bands), where conduction electrons have more 
states to be scattered into. States in this energy range 
are more affected by impurity concentration changes be- 
cause | Ae d \^>\ Ae s |. The most striking feature, how- 
ever, is the very pronounced maximum that appears in 
A -1 near the top of the d band, for reasonably large 
positive values of Ae d . The maximum becomes bigger 
and broader with increasing impurity concentration, and 
corresponds to a considerably faster decaying of g with 
£, as illustrated in the inset of figure 1(a). It arises from 
an interesting combination of disorder, hybridization and 
screening effects in the nanowire. To understand its phys- 
ical nature we recall that hybridization mixes the s- and 
d-conducting channels, and d-states are more directly af- 
fected by disorder in transition metal alloys. Screening 
in such systems is also very effective, and can cause sig- 
nificant changes in the local density of states (LDOS) 
around the impurity sites. It is well known that reso- 
nances and virtual bound states may appear near Ep, 
depending on the relative strength of the impurity po- 
tential. At sufficiently low temperatures, the conduc- 
tion electrons are scattered predominantly into available 
states near the scattering centers. Therefore, impurity 
resonant states at Ep may strongly influence electrical 
conductance, especially in disordered metallic systems. 
We notice that the energy position and some features of 
the A -1 maximum are closely related to the impurity po- 
tential strength. This is illustrated in the inset of figure 
1(b), where we compare the LDOS on a single impurity 
site, and its averaged value over the impurity neighbor- 
ing sites, with A ; all calculated in the same energy 
range for Ae d = 0.6, but smaller values of x = 0.05, and 
7 = 0.2, in order to highlight the virtual bound state 
structure in the LDOS. The peaks in A -1 and in the 
LDOS are clearly correlated in this case, as expected. 



0.08 




FIG. 1: Inverse localization lengths calculated for a non- 
magnetic disordered wire as functions of energy for different 
impurities (Ae d = 0.9 (a); Ae d = 0.6 (b)), and concentrations 
(x = 0.1 (dashed line); x — 0.2 (solid line)). The inset in (a) 
shows the average conductance (calculated for Ae d = 0.9 and 
x — 0.2) plotted against the wire length I for distinct values of 
energy: Ef = 2.4 (dashed line), Ef = 1.2 (dot-dashed line), 
and Ef = 0.8 (solid line). Lengths and energies are measured 
in units of the lattice spacing and s-hopping integral, respec- 
tively. The inset in (b) shows the LDOS (in arbitrary units) 
on an impurity site (dashed line), and averaged over its sur- 
rounding sites (thin solid line), together with A -1 (thick solid 
line), all calculated in the same energy range for Ae d = 0.6, 
but smaller values of x = 0.05, and 7 = 0.2. 



In fact, screening and disorder effects in ferromagnetic 
transition metal alloys can be highly spin dependent. In 
some Ni based alloys, screening is done largely by down- 
spin electrons, whereas up-spins normally play a signifi- 
cant part in Fe based alloys. Nii_a;Ti x , for instance, is a 
classical example where the introduction of a Ti impurity 
leads to an up-spin bound state being pushed completely 
above Ep JO)]]. The characters of up- and down-spin 
conducting states can certainly be altered by varying the 
alloy constituents and composition. To illustrate how 
appropriate combinations of elements may lead to very 
large spin filtering effects, we finally consider a ferromag- 
netic disordered nanowire attached to two non-magnetic 
metallic leads. We assume a typical value of exchange 
splitting Ajj = 0.5 on the d orbitals of the magnetic host 
atoms, and choose Ep — 1.2, corresponding to a strong 
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FIG. 2: Average conductances of majority- (solid line) and 
minority-spin electrons (dashed line), calculated as functions 
of wire length £, for ferromagnetic disordered nanowires. The 
calculations were performed for A d = 0.9, E F = 1.2 and 
x = 0.2 (see text). Lengths and energies are measured in units 
of the lattice spacing and s-hopping integral, respectively. The 
inset shows the percentage fraction of up-spin conductance 
P = gi/g plotted against I for Ef = 1.2 (solid line) and 
Ef = 0.35 (dot-dashed line), corresponding to a strong and 
weak ferromagnet, respectively 



ferromagnet with its up-spin d band completely filled. 
Our results for g^ and g±, calculated as functions of £, 
are shown in figure 2. They depict §i decreasing much 
faster than g^ , showing that the localization length of up- 
spin electrons in this case is much larger than that of the 
down-spins. Such a remarkable behavior is basically due 
to the large s character of the up-spin conduction states, 
contrasting with the much less mobile and predominant 
d character of the down-spin ones. As a consequence, 
the percentage fraction of up-spin conductance increases 
rapidly with the wire length. This is represented in the 
inset of figure 2, where P = g^/g is plotted against I for 
two different values of Ep , corresponding to a strong and 
weak ferromagnet, respectively. One immediately sees 
that very high spin polarized electrical currents can be 
achieved with relatively short ferromagnetic disordered 
nanowires, just a few hundred atoms long. This makes 
them excellent candidates for being used as spin injectors 
in nanoelectronic devices. 

In summary we have shown that localization lengths 
of up- and down-spin electrons may be rather differ- 
ent in ferromagnetic disordered nanowires. As a result, 
electrical currents flowing through such structures may 
acquire a high degree of spin polarization. The local- 
ization lengths for up and down spins are determined 
mainly by the character of the electronic states around 
the Fermi energy. By judiciously choosing the nanowirc 
alloy constituents and concentration, one can certainly 
modify such characters, and perhaps even control those 
localization lengths. There is a large ground for exper- 
imenting with such systems. Presently, the nanowires 



reported in reference M are composites of a carbon nan- 
otube in the core, covered by a relatively small amount 
of Ti, and a second metal on the outside. Nevertheless, 
careful heat treatments, and simultaneous deposition of 
different metals, may possibly be employed to fabricate 
nanowires of ferromagnetic disordered alloys. The most 
obvious choices to start with would be NiTi and FeTi 
alloys. The Ti impurity potential strength is relatively 
stronger in Ni alloys, but the magnetization of FeTi re- 
mains finite for larger concentrations of Ti. Other tran- 
sition metal combinations, may deserve to be examined. 
The fact that both Ni and Fe separately form continuous 
nanowires, suggests that NiFe alloys with virtually any 
concentration may also do, thus broadening the possibili- 
ties for exploring different spin dependent characteristics 
of the conducting states. NiCr alloys, with the Fermi 
level sitting on an up-spin virtual bound state, seems 
an interesting system too. The spin filtering mechanism 
reported here relies on localization and hybridization ef- 
fects. Most importantly, on the existence of relatively 
large contrast between the conduction state characters 
of up and down spins in the nanowires. We hope our 
findings will stimulate further investigation on these sys- 
tems. 
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